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We discuss the possibility of realizing metal-insulator transitions with ultracold atoms in two-dimensional 
optical lattices in the presence of artificial gauge potentials. Such transitions have been extensively studied 
for magnetic fields corresponding to Abelian gauges; they occur when the magnetic flux penetrating the lat- 
tice plaquette is an irrational multiple of the magnetic flux quantum. Here we present the first study of these 
transitions for non-Abelian 17(2) gauge fields, which can be realized with atoms with two pairs of degenerate 
internal states. In contrast to the Abelian case, the spectrum and localization transition in the non-Abelian case 
is strongly influenced by atomic momenta. In addition to determining the localization boundary, the momentum 
fragments the spectrum and the minimum energy viewed as a function of momentum exhibits a step structure. 
Other key characteristics of the non-Abelian case include the absence of localization for certain states and satel- 
lite fringes around the Bragg peaks in the momentum distribution and an interesting possibility that the transition 
can be tuned by the atomic momenta. 
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Ultracold atom research presents many avenues to explore 
problems at the forefront of theoretical physics. One recent 
line of investigation is the generation of artificial "magnetic" 
fields in an electrically neutral quantum gas JjJ which 
have the same effect on neutral atoms as a magnetic field has 
on charged particles. For atoms confined in an optical lat- 
tice, these synthetic gauge fields can be realized by impos- 
ing appropriate complex phase factors upon the amplitudes 
for tunneling of atoms between neighboring sites. This has 
renewed interest in the substantial theoretical literature of dia- 
magnetism of electrons in incommensurate two-dimensional 
lattices ||4j], a system which exhibits rich phase structure, but 
has seemed difficult to realize with any clarity in a condensed- 
matter system. Atoms in optical lattices may now provide a 
vehicle for experimental exploration of this subject. 

Furthermore, it also seems possible to create non-Abelian 
gauge fields in optical lattice systems |5], fields in which the 
components of the vector potential are noncommuting opera- 
tors rather than scalar functions. Such fields are of fundamen- 
tal importance in condensed matter and particle physics |6], 
and ultracold atomic physics now offers the prospect of con- 
trolled, well-characterized experiments on them @|: for ex- 
ample, the observation of analogue "magnetic monopoles" [8] 
and non-Abelian interferometry fl]. Recently, it was sug- 
gested that the artificial magnetic fields may provide a labora- 
tory realization of the fractal energy spectrum of an electron 
gas in a two-dimensional lattice, the Hofstadter butterfly |9], 
and the "Hofstadter Moth", the corresponding spectrum for 
the U(2) gauge fields ll|,|5[]. Another interesting feature that 
can be realized in these systems is the metal-insulator transi- 
tion when the magnetic flux a per lattice plaquette is an irra- 
tional multiple of the magnetic flux quantum [4]. In contrast 
to the Mott transition, which has its origin in particle inter- 
actions I Kill , this phase transitions is induced by competing 
periodicities associated with the cyclotron frequency and the 
frequency of motion of a particle in a periodic lattice [9]. 

The localization transition obtained by tuning the ratio of 



the tunneling along the two directions of the lattice has been 
extensively studied |4]. Here, we revisit the transition for the 
case when the "magnetic field" corresponds to U(2) gauge 
fields. We describe various novelties associated with the tran- 
sition in the non-Abelian case and discuss the experimental 
feasibility of seeing these transitions in two-dimensional opti- 
cal lattices. To best of our knowledge, this is the first study to 
examine the role of gauge groups on a quantum phase transi- 
tion. 

Our starting point is a tight binding model (TBM) of a par- 
ticle moving on a two-dimensional rectangular lattice (x, y), 
with lattice constants (a, b), and nearest-neighbor hopping 
characterized by the tunneling amplitudes (J, J A). When a 
weak external vector potential, A(x, y) = (A x , A y , 0), is ap- 
plied to the system, the Hamiltonian is given by 




where p is the momentum operator. For simplicity, we work 
in the Landau gauge: A(x, y) — (A x , A y (x), 0), where A x is 
a constant and A y depends only on the x-coordinate. The re- 
sulting Hamiltonian is cyclic in y and thus the eigenfunctions 
are proportional to plane waves in the y-direction. Denoting 
the transverse wave number of the plane wave as k y = k y /b, 
the wave function can be written as: ^(ma,nb) = e tkyn g m 
with x — ma and y = nb. For the Abelian gauge case, with 
A = (0, By, 0), one obtains the Harper equation 

9m+i + 9m-i + 2Acos(27ram - k y )g m = E g m , (1) 

where a = Bab/ (he/ e) and E(a) is the energy in units of 
J. The energy spectrum is the union over k y of individual 
energy spectra of the Harper equation. For rational a = p/q, 
the spectrum consists of q bands which are usually separated 
by gaps. As k y varies, the bands shift and their length may 
change, but they do not overlap, except at the band edges. For 
irrational a, the spectrum is independent of k y and the sys- 
tem exhibits a metal-insulator transition at A = 1. As we 
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discuss below, this insensitivity of the spectrum and the local- 
ization transition to the transverse wave number is lost when 
the gauge potential is non-Abelian. 

We obtain a non-Abelian vector potential for ultracold 
atoms in quite a natural way, when we treat systems in which 
the atoms have two degenerate internal states. The most gen- 
eral vector potential will couple the two states, and thus gives 
rise to U(2) gauge symmetry. In our treatment of the non- 
Abelian case, we follow the convention of an earlier study J2] 
and adopt its standard form of the vector potential, 

ea\2\ 1 -1 J ' V 2nma 2 )> )' 

Here ctfj determines the "magnetic flux" of the lattice. 

The wave function ^>(m, n) = e lkyn g m is a two- 
component vector g. m = [ ^ m J that satisfies the following 
equations, 

(t::) + Ci::) + (c.f)C:)-'^ 

where e£* = E — 2A cos(27rai7n — k y ). 

These coupled equations can be decoupled resulting in a 
pair of TBMs, 

+ = C m (a.\ , Ct 2 )9 m , (3a) 

e m-i»?m+2 + e m+i r im-2 = C m {a 2l a 1 )r\ m , (3b) 

where C m {a x , a 2 ) = CC+1C-1 ~ C-i ~ C+i- 

For ol\ = Oi2, the non-Abelian problem reduces to the 
Abelian problem described by the Harper equation [Eq. JJJ]. 
Calculation of the spectral properties and the detailed explo- 
ration of their localization transition in the non-Abelian prob- 
lem is done using various methods which include a transfer 
matrix approach 01, mapping to a dynamical system lHUl 
and the scaling arguments from Fibonacci renormalization 
theory 1 12]. In our study we fix ot\ to be the golden mean, 
a\ = (v5 — l)/2; in the transfer matrix analysis, a\ is ap- 
proximated by a sequence of Fibonacci numbers F n : F\ = 1, 
F 2 = 1, F n +\ = F n + F n -\. The non-Abelian problem is 
studied for various values of a 2 and we describe our results 
for a 2 — \,\ and a\. 

For the uncoupled Eq. ( I3al ( l3b» . the transfer matrix method 
is used to obtain the allowed energies of the system 00]. In 
general, Eqs. d3al(l3b» give two sets of allowed energies de- 
pending upon whether m is even or odd and therefore, the full 
spectrum of the non-Abelian system described by Eq. (0 is 
obtained by the intersection of the two spectra obtained for 
the even and odd-site 6 (or rf) equation. However, when otx 
and a 2 are both irrationals, or if one of them is zero, the spec- 
tra of the even and odd-site equations are identical. 

Localization properties, particularly the phase diagram 
of the system can be best studied by mapping the TBMs 
(Eq. d3al ) to a dynamical system, a discrete driven map ob- 
tained by defining x m = %— ^ 1 131 . The 9 equation becomes 



x m+2 = -l/[R m x m - f m ]. Here R m = C+i/ e m-i and 
fm = C m (ai, a 2 )/e" l 2 _ 1 . This mapping of a TBM to a dy- 
namical map is extremely useful in determining the localiza- 
tion transition as the absolute value of the Lyapunov exponent 
7 is related to the inverse localization length £ of the TBM, 
= 7/2. It should be noted that the maps are non-chaotic 
and hence the Lyapunov exponents are always negative. 
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FIG. 1 : (color online) Energy spectrum viewed as a function of k y 
for ai = T?7, 012 = 5 and A = 0.3. Light dots show the corre- 
sponding Abelian case for a± = a?2 = 89/144 where spectrum does 
not depend upon ky. Minimum energy shown by a solid line shows 

k 1 

the step structure with ^ = j being the global minimum or the 
ground state. 

In contrast to the Abelian case, the spectrum and the lo- 
calization properties of states in the non-Abelian case are 
strongly dependent on the transverse momentum k y . Figure^ 
shows the fragmented energy spectrum of the a 2 = | case as 
the transverse wave vector is varied. The fc y /27r = j emerges 
as a special wave number which characterizes the ground state 
and is also important for the central state (E = 0) as three dif- 
ferent bands touch at E = 0, ^ = |.[l3] Furthermore, the 
minimum energy develops a step structure and the number of 
steps increases with A (Fig. [3 and this structure resembles a 
fractal set at the onset to localization. The origin of these dis- 
continuities is two-fold: firstly, only a fraction of the full range 
of k y results in the allowed energies and secondly, the incom- 
patibility of the two sets of energy solutions corresponding to 
the even and the odd site equations. However, there exists a 
characteristic value of the wave number for which even-odd 
site equations result in identical set of allowed energies. For 
a 2 = i, this selected set of wave numbers is ^ = j (mod 
1). For a 2 = \, which exhibits analogous step-structure, the 
characteristic wave numbers are given by ^ = i (mod i). 

Figure [2] shows the localization aspect of the spectrum. In 
stark contrast to the Harper equation where all states localize 
simultaneously at A = 1, localization in the non-Abelian case 
begins at the band edges at A ps .48 for a 2 = \. In this case, 
all states except those in the central band localize before or at 
A = 1 . As A increases, the width of the central band decreases 
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and as A — > oo, only E(k y — it/ 2) — remains extended, a 
fact that can be inferred from Eq. (0 since it corresponds to 
e m — 0- This defiance of localization of the central band is a 
unique feature of the a 2 = \ case. For a 2 = \, localization 
begins at A w 0.24 and all states localize before A < 0.83. 
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FIG. 2: (color online) Upper panel shows minimum energy for 
a.2 = | with A = 0.4. Short and long dashed lines show the cor- 
responding Abelian case with a = ^ and |, respectively. Lower 

panel illustrates the localization of the ground state = |); lo- 
calization length f (in units of a) that with A = 0.4 (solid line) is 
oo becomes finite at A = 0.6 (dashed line). Here E is treated as a 
parameter: the extended states correspond to l/£ — > 0, spikes de- 
scribe localized states and the smooth lobes correspond to forbidden 
regions. 

To understand the origin of the selected set of wave num- 
bers, we now focus on the TBM corresponding to the E = 
state, which is always an eigenvalue of Eq. (|2jl for irrational 
a%. For a 2 = §, the & equation for E = reduces to, 

9 m+ 2 + m -2 + 2A e ffCOs(2%a 1 m-k y )6 m = (-2)6 m , (4) 

where A e f / = 2 A 2 cos(k y ). Thus, the E = state of the non- 
Abelian system is mapped to the E = — 2 state of a Harper- 
like equation @ where the strength of the onsite quasiperiodic 
potential is a sinusoidal function of k y . Diagonalization of 
Eq. @ for various values of A e ff shows that E — —2 is an 
eigenvalue of the system provided A e // < 0.48. Thus the 
boundary curve E(k y , A) — for the delocalized phase of the 
E = state is given by 2A 2 cos(k y ) w 0.48. This result is 
in perfect agreement with the numerically obtained boundary 
curve in Fig. fusing dynamical mapping and renormalization 
methods. 

A very direct evidence of the importance of the transverse 
wave number in the localization transition emerges for a 2 — 
j. The 9 equation for the E = state can be mapped to a 
non-hermitian TBM, 

Om+2 + 8m-2 + 2iA e cos(27raim - k y )6 m = 0, (5) 

where A e = 2A 2 cos(k y ) if m is odd and A e = 2A 2 sin(fc y ) 
if m is even and 9 m — i m 9 m . Analogous to the corre- 
sponding hermitian problem, the system exhibits self-duality 
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FIG. 3: (color online) Shaded region shows the extended phase for 
the E = state for 02 = a\. Long and short lines show the bound- 
ary of the extended phase for a.2 = \ and a.i = \, respectively. 
Solid vertical line shows the localization boundary for the Abelian 
problem described by Eq. Q with a=golden mean. 



at A e = 1 and this self-dual point describes the onset of 
localization 11411 . Therefore, the boundary curve describing 
the localization transition for the E = state corresponds 
to the minimum value of A consistent with the solutions of 
2A 2 cos(fc y ) = 1 and 2A 2 sin(A: y ) = 1. It should be noted 
that for ky/2ir = | (mod i), even-odd boundary curves are 
identical and we see the emergence of a characteristic wave 
number that corresponds to the maximum value of A for lo- 
calization of E = state. These analytic predictions are in 
perfect agreement with the numerical computation of the lo- 
calization boundary (See Fig.|5Jl. 

The above mapping elucidates the importance of the wave 
number k y in localization transition in the non-Abelian case. 
The transition threshold is a sinusoidal function of k y and 
there exists a special set of wave numbers that survive local- 
ization the longest as we increase the parameter A. The im- 
portance of k y and the existence of characteristic wave num- 
bers illustrated explicitly for E = state is also valid for 
ground states as it is precisely at these wave numbers that the 
full energy spectrum for even and odd sites for 9 (or 77) are 
fully compatible. This special set of wave numbers also corre- 
sponds to the ground state of the system. Interestingly, these 
features also characterize the localization transition for irra- 
tional a.2. The localization boundary for a 2 = ol\ and E = 
(See Fig. [3) clearly shows the emergence of a selected set of 
wave numbers k y /2n — j (n = 0, 1,2...) that remain ex- 
tended for the maximum range of A. Our detailed analysis of 
this case shows that the localization begins at A w 0.14 and 
all states are localized before A w 0.7. 

We now turn our attention to the experimental aspects 
of creating and observing the effects of non-Abelian poten- 
tials. A proposed experimental implementation of artificial 
Abelian |]l| and non-Abelian fields |5| consists of a two- 
dimensional optical lattice populated with cold atoms that oc- 
cupy two hyperfine states. The typical kinetic energy tun- 
neling along the s-direction is suppressed by accelerating the 
system or applying an inhomogeneous electric field in that di- 
rection. Tunneling is accomplished instead with laser driven 
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Raman transitions, induced by additional running wave lasers 
detuned to cancel the effect of the acceleration. To generate 
non-Abelian fields, we have two sets of Raman transitions. 
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FIG. 4: (color online) Momentum distribution for the ground state 
(ky/2-K — j) before it localizes for ai = and a-x = \ (A = 
0.4). The line with crosses shows the corresponding result for the 
Abelian case. Note that the central peak is truncated by a factor of 
three to show smaller peaks. 

We can experimentally access various A values (ratio of 
the tunneling in y and x directions) by adjusting the lattice 
beam intensity. Explicit formulas for tunneling due to kinetic 
energy (J) |15] and numerical calculations of laser-induced 
tunneling have been reported |1]. Based on those results, 
one can write A/_Er = (1/Ml) /(Vq/Eji, a), where / is 
a numerically known function. Here V$ is the optical po- 
tential of the lattice, En is the lattice recoil energy, and tl 
is the Rabi frequency associated with laser-induced tunnel- 
ing. In the non-Abelian case, there are generally two possi- 
ble values of A corresponding to tlx and tl 2 . By adjusting 
tl 2 /tli = f{V /E R ,ai)/f(Vo/E R , a 2 ), we obtain a single 
A in correspondence with the theoretical studies described 
here. For a practical range of Vq/Er between 5 and 45, 
/ decreases monotonically from 0.5 to nearly zero. To cre- 
ate a useful parameter range to see metal-insulator transition 
(0 < A < 2), we require that the parameter Ml /En be set to 
order unity. It is possible to achieve this with reasonable ex- 
perimental settings for alkali atoms like 87 Rb if the two states 
are taken to be the hyperfine levels of the 5 2 Si/ 2 state. To 
begin we note that the above scheme requires the Rabi fre- 
quency of the Raman transition tl, the detuning A, and the 
lattice trapping frequency v x to have well separated magni- 
tudes such that tl <C A -C v x to ensure that only the lowest 
band of the lattice is occupied and no other excitations oc- 
cur. Typical values of v x — En/h are on the order of tens 
of kilohertz HO]. The Raman transition is stimulated by two 
lasers with Rabi frequencies tlg 1 ^ and tli 1 ^ and a large de- 
tuning 5 r such that the effective Rabi frequency magnitude is 
tl = tig fie /2S r . The parameter Ml/ En can be fixed near 
unity if the Raman laser beams have intensities on the order 
of 1 mW/cm 2 (the intensity of one beam must be more than 
about ten times the other to satisfy the tl -C A condition) and 



are detuned from the 87 Rb D 2 line by about S r w 500 GHz. 

In cold atom optics, metal-insulator transitions can be de- 
tected by measuring the momentum distribution of the ex- 
panding atomic clouds. The Bragg peaks in the metallic 
phase (as shown in Fig. |4), reflecting the quasiperiodic or- 
der, become flattened at the onset to localization. It should be 
noted that in the non-Abelian case, the momentum distribution 
shows additional satellite peaks around the Bloch vectors. 

In this paper, we have identified means for driving a lo- 
calization transition in two-dimensional optical lattices, and 
have contrasted the effects of different gauge symmetries. The 
U(2) gauge system clearly offers a much richer environment 
than does £7(1). For example, using polarized fermions such 
as 40 K and with the Fermi energy in the central band, one 
may observe the complete defiance of localization. For a BEC 
loaded in the lattice, smooth changes in the ground state en- 
ergy for the Abelian system will be replaced by extremely 
sharp oscillations as the transverse component of the momen- 
tum is tuned (see Fig.[0. In addition to controlling the nature 
of the quantum phases by the ratio of the tunneling in two di- 
rections, the transport properties can also be controlled by the 
momentum. Such novelties will broaden the possibilities of 
using ultracold gases in exploring new horizons in physics. 
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